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1 Motivations 


3 Computations of A s and A* 


The moist-air entropy is defined in Marquet (2011, 
hereafter Mil) by s = s re f + c p d ln(0 s ) in terms of 


two constant values (s re y, c p( f) and a potential entropy 
temperature denoted by 6 S . It is shown in Mil that 
a quantity denoted by (0 s )i plays the role of a leading 
order approximation of 6 S . 


The aim of this note is to demonstrate in a more 
rigorous way that (0 s )i is indeed the leading order ap¬ 
proximation of 0 S , and to derive a second order ap¬ 
proximation which may be used in computations of 
values, gradients or turbulent fluxes of moist-air en¬ 
tropy. Some impacts of this second order approxima¬ 
tion are described in this brief version of a note to be 
submitted to the QJRMS. 


2 Definition of 9 S and (+)i 


Let us define A s by the formula 6 S — 0i exp(A s q t ) 
where 0 S , 0i and q t are known quantities, and thus by 

As = ^ ln (+)' (3) 
The aim is to compute A s from © for a series of 16 

vertical profiles of stratocumulus and cumulus with 
varying values of qt , 6 S and 0/, in order to analyse the 
discrepancy of A s from the constant value A r ~ 5.87 
given by ©• 

A trial and error process has shown that plotting 
A s against In (7+) leads to relevant results (see Fig{l]). 
Clearly, all stratocumulus and cumulus profiles are 
nearly aligned along a straight line with a slope of 
about —0.46. It is likely that this slope must corre¬ 
spond to — 7. This linear law appears to be valid for 
a large range of r v (from 0.2 to 24 g kg -1 ). 


The potential temperature 6 S is defined in Mil by 


0 S = {0 s) 1 





(I+7777) 
(1+7777 .) KSqt 


where 


{0s) 1 = 0i exp(A r q t ) . 


(1) 


This definition of 0 S is rather complex, but the 
more simple quantity (0 s )i was considered in Mil 
as a leading order approximation of 0 S , where the 
Betts’ potential temperature is written as 9i = 
6 exp[ - (L v qi + L s qi) / (c pd T )]. 


The total water specific content is qt = q v + qi + 
qi and r v is the water vapour mixing ratio. Other 
thermodynamic constants are: Rd ~ 287 J K -1 , R v ~ 
461.5 J K -1 , c pd « 1005 J K -1 , c pv « 1846 J K -1 , 
k Rd/c p d ~ 0.286, A c pv / c pd - 1 « 0.838, 5 
R v /Rd—1 ~ 0.608,77 = R v /Rd ~ 1.608, , e — Rd/Rv ~ 
0.622 and 7 = R v /c p d ~ 0.46. 


The term 

A r [ ( S v ) r ( £>d)r ] / Cpd ~ 5.87 


( 2 ) 


depends on reference entropies of dry air and water 
vapour at T r — 0 C, denoted by (sd) r — s d(Tr 5 c r ) and 
(5^) r = s v (T r ,p r — e r ), where e r = 6.11 hPa is the 
saturating pressure at T r and p r — 1000 hPa. The two 
reference entropies (s^) r ~ 12673 J K -1 and (sd) r ~ 
6777 J K -1 are computed in Mil from the Third Law 
of thermodynamics, leading to A r ~ 5.87. The refer¬ 
ence mixing ratio is r r — ee r / (p r — e r ) ~ 3.82 g kg -1 . 



Figure 1: A plot of A s against ln(r v ) for 8 cumulus (dashed 
red), 7 stratocumulus (solid blue) and ASTEX (solid black) 
vertical profiles. The constant value A r ~ 5.87 corresponds 
to the horizontal dashed black line. An arbitrary line with a 
slope of —0.46 is added in dashed-dotted thick black line 

It is then useful to find a mixing ratio r* for which 
A s « A* = 5.87 - 0.46 In (r„/r*) (4) 

hold true, where r* will play the role of position¬ 
ing the dashed-dotted thick black line of slope — 7 « 
—0.46 in between the cumulus and stratocumulus pro¬ 
files. This corresponds to a linear fitting of r v against 
exp[ (A r — A s )/7 ], r* being the average slope of the 
scattered data points. It is shown in Figj2] that the 
value r* « 12.4 g kg -1 corresponds to a relevant fit¬ 
ting of all cumulus and stratocumulus vertical profiles 
for a range of r v up to 24 g kg -1 . 





















Then, the first order expansion of (l+rjr v )^ K A+ S( it) 1 
for small 77 ~ qt is equal to exp( 7 q t ), since 7 = tc 77 . 
The last term (1 + 77 r r )( K s qt ^ leads to the higher order 
term exp[ 7 5 q t r r ] ~ exp[ O(q^) ], since r r < 1 and 
gt « 1. Other terms depending on temperature and 
pressure are exactly equal to exp[ A q t In (T/T r ) ] and 
exp[ - k, 5 q t In (p/p r ) ]. 

The Taylor expansion of 9 S can thus be written as 


9 S « 6 exp | 


^vap Ql T Lsub qi 
c p d T 


exp (A* q t ) 


(5) 


Figure 2: Same as Fig\l\ but with r v plotted against the 
quantity exp[ (A r — A s )/7 ]. Two lines of slope r* = 10.4 
and 12.4 g kg -1 are added as dashed-dotted thin black lines. 


x exp < q t 


A In 


-kS Inf— 

\PrJ J 


+ o(« 


where 


A* = A r - 7 ^( 77 / 77 ) 


and 77 = t> x exp(l) ~ 10.4 g kg 1 (see Figs{2]and[3]). 


It is shown in Figj3]that A s can indeed be approx¬ 
imated by A*( 77 , 77 ) given by Q, with a clear im¬ 
proved accuracy in comparison with the constant value 
A r ~ 5.87 for a range of r v between 0.2 and 24 g kg -1 . 
Curves of A*( 77 , 77 ) with 77 = 10.4 and 12.4 g kg -1 
(solid black lines) both simulate with a good accuracy 
the non-linear variation of A s with 77 , and both simu¬ 
late the rapid increase of A s for 77 < 2 g kg 1 . 



Figure 3: Same as but with r v plotted against the 

quantity A s . The two thin black lines correspond to 0 with 
77 = 10.4 or 12.4 g kg ~ x . 


4 Mathematical computation of A* 

It is important to confirm, by using mathematical ar¬ 
guments, that ( 0 S ) 1 corresponds to the leading order 
approximation of 0 S , and that the slope of — 7 « —0.46 
(with 77 ~ 10.4 or 12.4 g kg -1 ) corresponds to a rele¬ 
vant second order approximation for 0 S . These results 
are briefly mentioned in Marquet and Geleyn (2015). 

First- and second-order approximations of 6 S can be 
derived by computing Taylor expansions for most of 
terms in the first formula recalled in Section [2l 

The main result is that the term (r r /r v )^ qt ^ is ex¬ 
actly equal to exp[ — (7 q t ) In( 77 / 77 ) ]. 


The first order approximation of 0 S is thus given 
by the first line of ^ with A* = A r , namely by the 
expected (0 s )i. An improved second order approxima¬ 
tion is obtained by using A* instead of A r and by taking 
into account the small corrective term — 7 ^( 77 / 77 ). 

Impacts of the second line of © with terms de¬ 
pending on temperature and pressure lead to higher 
order terms which must explain the fitted value 77 ~ 
12.4 g kg -1 observed for usual atmospheric conditions. 

5 Conclusions 

It has been shown that it is possible to justify math¬ 
ematically the first order approximation of 6 S derived 
in Mil and denoted by (0 s )i ~ 0i exp (A r q t ), which 
depends on the two Betts’ variables (0i,qt) and A r « 
5.87, leading to s« s ref + c pd In(0 ; ) + c pd A r q t . 

A second order approximation is derived and com¬ 
pared to observed vertical profiles of cumulus and stra- 
tocumulus, leading to s ^ s re f + c p y ln(0j) + c p g A* q t , 
where the second order term A* = A r — 7 In( 77 / 77 ) 
depends on the constant A r ~ 5.87, on the mixing 
ratio 77 , and on a tuning parameter 77 ~ 12.4 g kg -1 . 

The use of the second order term A* depending on 
the non-conservative variable 77 ~ q t — qi — qi can 
explain why it is needed to replace the Betts’ potential 
temperature 61 for computing flux of moist-air entropy: 

w'0' s « exp (A* q t ) w'0[ + A* 0 S w'q' t - (7 q t 0 s /r v ) w'r ' v , 
the last term depending on w'q' t minus w'q[ or w'q '. 
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